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Abstract 

^ We obtain analytic formulae for the null geodesies of Friedmann-Lemaitre-Robertson- 

Walker spacetimes with scalar perturbations in the longitudinal gauge. From these we 

> 

^ ' provide a rigorous derivation of the cosmological lens equation, and obtain an expression for 
O ' the magnification of a bundle of light rays without restriction to static or thin lens scenarios. 
O , We show how the usual magnification matrix naturally emerges in the appropriate limits. 

o 

^ • 1. Introduction 
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^ ! The bending of light by a single symmetric gravitational lens in Euclidean space is 

\ shown in Fig. 1. The symmetry guarantees that the lines of sight from the observer, o, 
to the lens, /, and to both the lensed and unlensed images of the emitter, e, are coplanar, 
and can be related by the angles a, (3, and 9. D{oJ) and D{o,e) are the distances from 
the observer to the lens and emitter, respectively. We assume that the defiection angle, a, 
is small. Then, locally about the line of sight to the emitter's image, we may approximate 
the two-spheres at distances D{oJ) and D{o,e) from the observer as planes, called the 
lens plane and source plane respectively. D{1, e) denotes the distance between these two 
planes. The assumption of small defiection angle also allows us to relate the lensing angles 
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by P = 9 — aD{l,e)/D{o,e). This is the simplest expression of the gravitational lens 
equation (Refsdal 1964). 

The generalization of this equation to more complicated lens structures and non- 
Euclidean background spaces proceeds by a number of steps. We continue to assume 
that the actual path of the photon is well approximated by two segments joined at a point 
of deflection, p, located near the lens. A general lens is not symmetric so that the angles a, 
P, and 9 are not necessarily coplanar. To handle this, consider a set of cartesian axes with 
origin at the observer. Choose the a;-axis to coincide with the line-of-sight to the image. 
Let a*, where i runs over {2, 3}, be the angle between the a;-axis and the projection into the 
a;*a;-plane of the line-of-sight vector from the deflection point, p, to the emitter. Similarly, 
let be the angle between the projection into the x*x-plane of the line-of-sight vector from 
the observer to the lens and the projection into the a;*a;-plane of the line-of-sight vector 
from the observer to the unlensed image. Also, let 9^ be the angle between the x-axis and 
the projection into the a;*a;-plane of the line-of-sight vector from the observer to the lens. 
To allow for non-Euclidean spatial geometries D{o, e) is taken to be the angular-diameter 
distance in the background geometry from the observer to the intersection of the x-axis 
with the source plane, and D{l,e) is taken to be the angular-diameter distance in the 
background geometry between the deflection point, p, and the intersection of the cc-axis 
with the source plane, p'. Then, again assuming small deflection angle. 

This is the standard cosmological gravitational lens equation (e.g. Schneider, Ehlers, and 
Falco 1993, Chapter 2). An important quantity associated with this formalism is the 
magniflcation matrix, M*^ = dP'^/d9^, which contains information on the deformation 
of ray bundles connecting the observer and emitter. For example, the inverse of the 
determinant of this matrix is the magnification of an image relative to an unlensed image. 
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The purpose of the current paper is to address a number of subtle issues that arise 
when we attempt to justify mathematically the use of the lens equation (1) for calculations 
in our Universe, although most workers agree that the physical basis for its use in observed 
lens systems is strong. First, there is the question of the correct choice of distance factors. 
There exists a large literature addressing this question, primarily concerned with the ap- 
propriateness of the so-called Dyer-Roeder distances (Dyer and Roeder 1972, 1973; Ehlers 
and Schneider 1986; Futamase and Sasaki 1989; Watanabe and Tomita 1990; Watanabe, 
Sasaki, and Tomita 1992; Sasaki 1993). Our work suggests that within the framework 
of cosmological perturbation theory, the natural distance factors to use are those of the 
background. Hence, the choice of distance factors is equivalent to the choice of cosmolog- 
ical model, in agreement with the recent results of Sasaki (1993). The issue of the most 
appropriate choice of cosmological model must be addressed in its own right. 

The second concern in any mathematical investigation into the lens equation is the 
accuracy of the approximation of an actual photon path by two geodesies of the background 
which join at a point near the lens: the deflection point, p. On physical grounds we 
expect this approximation to be good for systems for which the photon-lens interaction is 
localized: the thin-lens approximation. One purpose of our present work is to quantify the 
error involved in using two geodesies of the background, rather than the actual path, in 
deriving the lens equation, (1). 

Third, how are we to find the angles appearing in the lens equation from physical data? 
Generally, the a* are taken to be those calculated in Einstein-de Sitter spacetime, since the 
overall curvature of space should not be important near p, where the light ray interacts with 
the lensing object. For static, thin lenses, the deflection angle is written as a superposition 
of point mass deflection angles contributed by mass elements of the lens projected onto 
the lens plane (Schneider et al. 1993). For brevity we will term the resultant angle the 
"Einstein angle." Another purpose of the present work is to derive this result from the full 
equations of light propagation under an appropriate set of mathematical approximations. 
Succinctly, the lens equation, for static, thin lenses, effectively assumes that the light path 
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is described by the Jacobi equation of the background spacetime subject to an impulsive 
wavevector deflection at the lens plane by an angle equal to the usual Einstein bending 
angle. We wish to quantify the level of approximation involved. 

There have been two notable recent attempts to clarify the validity of the cosmological 
lens equation by deriving it from the optical scalar equations (Seitz, Schneider, and Ehlers 
1994) and the Jacobi equation (Sasaki 1993). However, a crucial difference between these 
papers and the present work is that they treat the path of the light ray differently near to 
and far from the lens. It is precisely this assumption that we must eliminate if we hope to 
gain a more general lens equation able to quantify the errors implicit in equation (1). 

Our approach is to investigate the cosmological lens equation as it emerges from the 
geodesic equation. In this our work is complementary to that of Seljak (1994) and Kaiser 

(1992) who have used an approach like this to investigate certain lensing questions for a per- 
turbed Einstein-de Sitter spacetime. Where our work overlaps that of these authors we are 
in agreement. We will show, however, that it is possible to handle the curved Friedmann- 
Robertson- Walker (FRW) spacetimes by analogous calculations, though of somewhat more 
technical difficulty. With theorists beginning to take the idea of spatially curved models 
more seriously (Kamionkowski et al. 1993, 1994; Spergel et al. 1993), we feel that this 
extension is of more than formal importance. In this paper we make use of a method for 
constructing null geodesies in perturbed spacetimes introduced in Pyne and Birkinshaw 

(1993) , hereafter PB. The method is the analog for geodesic curves of familiar perturba- 
tion techniques for differential equations. The results presented here come from applying 
this method to FRW spacetimes with scalar perturbations in the longitudinal gauge. Our 
principal results are: 

(1) analytic formulae (equations (13) and (14)) for light rays in the spacetime (5); 

(2) a general expression for the magnification undergone by a bundle of light rays capable 
of handling non-static, geometrically thick, lenses (equations (35), (36), and (38)); 

(3) a demonstration that the usual lens equation (1) and magnification matrix are recov- 
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ered in the appropriate approximations. 

To our knowledge, this is the first rigorous derivation of the cosmological lens equation for 
perturbed FRW spacetimes with spatial curvature. 

The outline of this paper is as follows. In section 2 we review the perturbative geodesic 
expansion introduced in PB. In section 3 we apply this formalism to the problem of con- 
structing null geodesies of FRW spacetimes with scalar perturbations. With the help of the 
equation of geodesic deviation for the FRW background, the solutions we find will enable 
us to understand the role of the Einstein angle for light propagation in these spacetimes. 
In section 4 we obtain an expression for the magnification of a bundle of light rays in such 
spacetimes without restriction to static, geometrically thin perturbations. We then discuss 
the emergence of the usual cosmological lens equation, (1), in the appropriate limits. In 
section 5 we illustrate the use of the magnification equation by solving it for a point mass 
embedded in an Einstein-de Sitter spacetime. 

2. The Perturbative Geodesic Expansion 

We work in geometrized units, G = c = 1. We let Greek indices z^, ... run over 
{0, 1, 2, 3} and Roman indices ... run over {1, 2, 3}. The spacetime metric is taken to 
have signature +2. Our Riemann and Ricci tensor conventions are given by [Aq,, A^] = 

R^vapv" and Rap = R^a^i/S- 

In PB we showed that given a null geodesic of gj^J^ x^'^^^{X), with A an affine parameter, 
we could construct a set of four functions, which we call the separation, a;*^^^'^(A), trans- 
forming as a vector under infinitesimal co-ordinate change, which were sufficient to ensure 
that x^^^^^{X) +x^^^^^{X) is an affinely parametrized null geodesic of gj^J + h^i,. In order to 
construct x^^'>^{X) we need three quantities, two referring solely to the background met- 
ric and one referring to the perturbation. First, we need the parallel propagator (Synge, 
1960) appropriate to x^^^^{X) in gjj^J (in PB we called this the connector after DeFelice 
and Clarke (1990, Section 2.3), mostly to avoid confusion with the Jacobi propagator be- 
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low, but we prefer Synge's term). We remind the reader of one crucial property of the 
parallel propagator: if (A2) is a vector at x^^^^^ (A2) and (Ai) is its parallel trans- 
late along x^°)'*(A) to the point x^^^^ (Ai) then the parallel propagator P (Ai, A2)'^ i, obeys 
(Ai) = P (Ai, A2)'^ i/f ^ ('^2)- That is, it parallely propagates vectors along x^^\X). 

The second quantity we need, the Jacobi propagator, also refers only to the background 
spacetime. Introduced in PB, it is an 8 x 8 dimensional matrix which serves as a Green's 
function for the Jacobi equation for gj^J along a;'^'^(A). Its explicit construction makes use 
of the parallel propagator, the curvature tensor of gj^J, R^^^^i^pa, and the tangent vector 
to a;(°)'^(A), the wavevector, A;(°)'^(A) = dx^^^'^{X)/dX. With eight dimensional systems 
the usual tensor notation can be cumbersome so we will use the matrix notation of PB. 
We let 7^(A)^<^ denote the 4 x 4 matrix R^^'^^'^p^k^^^'^k^^'^P evaluated at x^^\X) and write 
P (Ai, A)'' {Xy pP (A, Ai)^ ^ as P (Ai, A) U (A) P (A, Ai). Then the Jacobi propagator 
U (Ai, A2) is given by 



u^>^-^>^-) = ^-^AC{p{XuX)n\x)p{XM o)^^)- 



Here 1^ is the 4x4 identity matrix and V denotes the path ordering symbol. 

The final quantity we need encodes the effects of the perturbation itself. It is a vector 
field defined along x^^^^{X). We denote the covariant derivative with respect to g^^^ by a 
semi-colon. Then the perturbation vector at affine parameter A, f^{X), is given by 



evaluated at a;^°)'^(A). 

It was shown in PB that the force vector, Jacobi propagator, and parallel propagator 
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allow us to solve for the separation, x^^^^(X). In particular, we suppose that we desire to 
construct a geodesic of gj^J + h^^ of the form x^{X) = x'^^^^^{\) + x^^^^{\) and that we 
know the appropriate boundary conditions at affine parameter Ai. Then x^^-*^ is given at 
arbitrary affine parameter, A2, by 



d 

d\2 



P(Ai,A2)x(i) (A2) 
P(Ai,A2)x« (A2) 



= C/(A2,Ai) 



_d_ 
dX 



x(i)(Ai) 
P(Ai,A)x«(A) 



+ 



J Ai 





i^(Ai,A)/(A) 



dX 



(4) 



where the integral in this equation is taken over the background path, x^^^^{X). 

While the solution above looks complicated, in practice background spacetimes are 
chosen specifically for their high degree of symmetry and tractability and this often allows 
us to construct the needed propagators explicitly. In these cases, equation (4) reduces the 
work of finding null geodesies in the perturbed spacetime to a simple problem of integration. 
We will see below that the crucial FRW spacetimes belong to this class. 

The consistency criteria for our solution are precisely those of all perturbative type 
geodesic calculations; extremely heuristically, the background and the constructed geodesic 
should not be allowed to reach regions where their spatial or temporal deviations are such 
that the geodesies effectively feel different gravities at equal affine parameter, either due 
to the perturbation or the curvature of the background itself. A simple a priori estimate 
for the domain of validity was given in PB, but in practice, the consistency may usually 
be checked easily after a solution is obtained. ^ 

^ We note that the condition <^ k, given in PB is not, in fact, necessary. We thank 
Uros Seljak for pointing this out to us. 
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3. The Deflection Angle 

Wc will now employ the techniques reviewed above to investigate the theory of grav- 
itational lenses in cosmology. Our starting point is a choice for the metric. We choose to 
work with FRW spacetimes with scalar perturbations written in the longitudinal gauge, 



where 7 = l+K.r'^/A, k being the spatial curvature parameter (±1 or 0) and = x^+y^+z^. 
Inhomogeneities are represented by the quasi-Newtonian potential, 0. For the order needed 
by us, the expansion factor, a(r/), is unperturbed from its Friedmann form (Jacobs, Linder, 
and Wagoner 1993). We choose this form for the metric for a number of reasons. The 
metric (5) is also used by Seitz et al. (1994) and Sasaki (1993), so that direct comparison 
of results is possible, and recent work by Futamase (1989) and Jacobs et al. (1993) has 
shown that structure of galactic scale and greater in our Universe can be well modeled by 
metrics of this type. While results obtained with the metric (5) are not appropriate for 
lensing by gravitational waves or vector perturbations, these cases are easily handled in 
an analogous manner. 

Next we choose a particularly useful class of background light rays with which to build 
our perturbed solutions; radial null geodesies intersecting the observer at the spatial origin. 
Because the Friedmann expansion, a, plays the role of a conformal factor it is simplest to 
work with the null geodesies of ds^, defined by ds^ = aP'ds^. Light rays in these two metrics 
coincide and their (affine) parameterizations are related by k'^ = a~'^k'^. With the observer 
located at the spatial origin, the radial null geodesies of ds*^*^^^, (i.e. of that part of ds'^ 
independent of 0), may be written /c^*^^^ = 1 and /c*-*^-** = —76*, where e* are the direction 
cosines at the observer, so that Z^f=i(e*)^ = 1 (McVittie 1964). Note that we have choosen 
our aflfine parameter to coincide with conformal time. This is purely for convenience. The 
explicit solutions for the comoving radius and for 7 along such rays are given by 





8 



r(A) = 2tan. (^V^) 
7(A) =sec^ 1 ^° ] 



where Aq is the affine parameter at the observer. The subscript k on a trigonometric 
function denotes a set of three functions: for «; = 1 the trigonometric function itself, for 
K = — 1 the corresponding hyperboUc function, and for k — the first term in the series 
expansion of the function. The paths of the rays are x^^^^ = A, x*^^^* = re*. 

The equations of parallel transport along x^^^^{X) are easy to solve. Two vectors, V2 
and Vi, at x^^"^^ (A2) and x^'^^^ (Ai) respectively, are related by parallel translation provided 
that V2 = Vi and that i'2 = 7 (A2) ^'1/7 (Ai). We can, thus, read off the parallel propagator 
for our class of geodesies. 



P{\2,Xir.= (l 7(A.').iV (7) 

7(Ai)''iy 



Next we will obtain the Jacobi propagator. The Riemann tensor of ds^^^'^ is non- 
vanishing only when all indices are spatial, when 



where gj^^J is the metric ds^^^'^. As a result we can write R^^^^i,pfjk^^^^k'^^^P = —kJ^^^j where 
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0, 

0* 5] - e'ej 



) 



(9) 



J is idempotent, = J. This allows us easily to sum the series defining the Jacobi 
propagator, (2), in 4 x 4 subblocks. The result is, returning fully to our matrix notation. 



We note that J may be interpreted as a projection operator into the space transverse to 
the photon direction in the comoving spatial hypersurfaces. Thus the Jacobi propagator, 
U, has split into a transverse rotation and a longitudinal shear. 

Another crucial quantity we will need before we can construct the perturbed light rays 
is the force vector f^^. We note that the force vector may be constructed not only from 
(3), but also from the equivalent 



either the metric perturbation or its first partial derivatives. Direct calculation yields 




(10) 



(11) 



where T^^^/^^fi denotes that part of the Christoffel connection of ds'^ which is linear in 




(12) 
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where a comma denotes an ordinary partial derivative. 

It remains only to determine the appropriate initial conditions for the separation before 
we can use (4) to gain the null geodesies of ds'^. We will choose the perturbed geodesic and 
the background geodesic to coincide at the observer, x^^^i^ (Aq) = 0. As discussed in PB, we 
can not take the wavevectors of the two geodesies to coincide fully at the observer because 
the geodesies must each be null in their respective, different, metrics. We will choose, for 
convenience, with k^^^^ — dx^^^^/dX, k^^^'^ (Ao) = 0. The constraint that k^^^'^ + k^^^'^ be 
null in ds"^ at the observer then tells us that k^^^'^ (Ao) = — 20o where we denote the value 
of at the observer by 0o. 

It is now simply a matter of assembling the necessary pieces in (4) to gain the separation 
at arbitrary affine parameter Ag, and thus the light rays of ds^. Some straightforward labor 
yields 



(Ae) = -2 (Ae - A^) 4>o + 2 dX (A - Ae) /c^^^^^^.m (A) (13) 



and 



(Ae) = -2/c(°)^ dX (A - Ae) |^(A) 



dr] 

+ 27(Ae) dAsin« (A - Ae) ^ {V±(Py (A) 
JXo 7 A 



(14) 



7(A) 



where V_l* = g(^)^'^ — e^e*^ dm is the transverse gradient operator. We see that the 
spatial separation is naturally written as the sum of a longitudinal and a transverse term. 
We remind the reader that the above integrals are taken over the background geodesic 
a;(0)M(^), Hence if these solutions are used for geometrically thick lenses, the error terms will 
become important at some finite distance along the geodesic. In this case, it will in general 
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be necessary either to apply an iterative procedure, incorporating a number of background 
paths, or to appeal to statistical arguments to bound the errors. These approaches are 
familiar from the usual multiple lens plane theory (Schneider et. al. 1993; see also Seljak, 
1994). A crucial difference, however, between the use of multiple background paths in our 
formalism and the multiple lens plane method is that formulae (13) and (14), in principal, 
allow the photon path to be approximated to arbitrary accuracy by succesive modifications 
of x^^^^{X), in contrast with the multiple lens plane theory where the continuum limit is 
not compatible with the assumptions underlying the theory. 

It is not too difficult to check by straightforward calculation that the geodesic defined 
by (13) and (14) is, in fact, null. This is to be expected in light of the general theorem 
proved in PB that our constructed geodesic preserves its null character. We point out, 
however, that our appeal, in that paper, to the co-ordinate invariance of scalar quantities 
in order to argue for the vanishing of the term involving the partial derivative of gjj^J is 
invalid. The theorem is, nevertheless, true. ^ 

We can gain more understanding of (13) and (14) by considering their relation to the 
equation of geodesic deviation. In Appendix A we show that the Jacobi equation of ds^^^'^ 
subject to an arbitrary impulsive wavevector perturbation Sk^ at some affine parameter 
u, is solved by deviation vector 5x^ with spatial components 



Sx'{X) = sin«(w - A) 6ki{u) - ^ {u - A) SkUu) (15) 

7(tt) " 

where Sk\^ = (Sj — e^ej^ 5k^ is the impulse in the transverse direction, and Skl^^ = Sk'^ — Sk\^ 
is the longitudinal impulse. 

^ One way to patch it up, for instance, is to note that gj^J-p = allows us to replace the 
offending term by a sum of two terms linear in the Christoffel symbols of the background, 
each of which vanishes in the chosen co-ordinate system of the second half of the proof. 
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A comparison of (14) and (15) leads to the interpretation of the spatial components of 
our solution for the separation, a;^-*^)*, as the result of a continuous sequence of impulsive 
perturbations 



Sk' = -2{V±(l)y + 2k^^^'^. (16) 



The form of the impulse can also be gained directly from our equations by differentiating 
the spatial separations with respect to the afhne parameter, and inserting a delta function 
at Xi into the integrand which forces the integrand to vanish except at the lens plane. This 
gives 



{Xi) = -2 (V^<^)^ {Xi) + 2A;(«)^ (A^) ^ (A,) , (17) 

which is exactly the impulse found above. 

At this point we need only do a little work to recover the Einstein angle from our 
equations. Establish a set of cartesian axes at the observer and choose the unperturbed 
wavevector k^^^ = (1,-7,0,0). Consider a static, localized perturbation in the xy-plane. 
Then the angle represented by the impulse perturbation found above (17), that is, the 
angle A;^*^)* (A;) + k^^^'^ (A;) makes with k^^^'^ (A^), is given by 



&^==^^^ = -2-,(X.)^„. (18) 

l{Xi) 



The factor of 7 (A^) is present only because our co-ordinates are scaled in an unusual way 
at the lens plane. If we make the co-ordinate change x''' ~ x'^/'y{Xi) the metric on the 
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lens plane becomes Minkowskian and, locally near the deflector, y' and z' serve as normal 
co-ordinates on the lens plane. In these co-ordinates (j)^yt on the lens plane takes on the 
usual Newtonian form (with origin shifted from the lens, accounting for the unusual minus 
sign). Since our lensing angle —2'^{\i)(^^y = —20^^/ the integrated impulse lensing angle 
for a localized perturbation is exactly the Einstein deflection angle. We emphasize that 
this is the flrst time this has been shown rigorously for the curved FRW spacetimes. 

For completeness we note that the timelike component of the separation can also 
be analyzed by comparison to the Jacobi equation. The Jacobi equation of ds^'^^'^ for an 
impulse wavevector perturbation Sk^ at affine parameter u results in a timelike component 
of the deviation vector 5x^ = —{u — X)Sk^{u). Comparison with our solution for the 
separation reveals that the time delay undergone by the light ray relative to the fiducial 
background ray may be considered to result from a sequence of impulses 5k^ = —2k^^^'^(f)^rai 
in addition to a boundary term. 



4. The Magnification 

We now want to examine the magnification undergone by a bundle of light rays. We 
define this after Schneider et al. (1993) in the following way. Suppose a source of given 

physical size at some redshift is observed to subtend solid angle dO. An identical source 
observed at identical redshift placed in an FRW spacetime would subtend solid angle dVL^'^^ 
The magnification M is defined to be dQ,/dVl^'^\ 

To gain the magnification we will construct an infinitesimal bundle of light rays in ds'^ 
which emanate from a source and converge at an observer, located at the spatial origin 
of co-ordinates, by varying the direction cosines of the background ray, e*, in equations 
(13) and (14). We will determine the solid angle, dO,^ subtended by the rays in the rest 
frame of an observer with four- velocity li^ = (l/tto) (l — 4'o,vlj- We will then ask what 
local area transverse to its direction of propagation the bundle sweeps out in a frame with 
four velocity tt^ = (l/og) (l — (j)e,vl), dA, at a given redshift, z. We regard the peculiar 
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velocities vl and vl as first order quantities so that and uf^ are properly normalized 
to first order in ds'^. We suppose a source with four- velocity uf^ intersects our bundle, 
that its redshift is z, that its physical size is dA, and that its shape such that it exactly 
fills the beam; that is, our bundle is the light of a physical source. In this way we will 
gain a relationship between dQ and the redshift and proper size of the source and the 
four-velocity of the observer. A similar relationship is easy to derive for an identical source 
in the background spacetime. Comparison of the two expressions will yield M. Figure 2 
illustrates the constructions of this section. 

We begin by choosing a set of null geodesies of the background with which to construct 
our congruence. We use the two-parameter family of curves given by x^^^'^ = ^A, re*) with 
r as in (6), = {l,dsm9,dcos9), d G (0, e) with e infinitesimal, and 9 G (0, 27r]. We work 
to first order in e. To this order, this set of rays defines a null congruence of ds^^^'^. To 
each of the rays of this congruence is associated a null ray of ds'^, and hence of cZs^, by (13) 
and (14) above, x^^ (A; d, 9). In the rest space of the observer, these rays define a cone. To 
see this, we note that the conical shape is clear for a comoving observer. For an observer 
with some peculiar velocity the assertion then follows from a result of Sachs (1961) on the 
geometry of null rays. 

Next we will determine the two-dimensional projected area of our bundle at a given 
redshift. In Appendix B it is proven that, to first order in the perturbation, = x^{d = 
e,9) — x^{d = 0) is a one-parameter family of geodesic deviation vectors of ds^ along 
the central geodesic of the bundle, x^^{d = 0). From this point on, any wavevector or 
path pertaining either to the perturbed or unperturbed spacetimes not written with an 
explicit (i-argument is intended to refer to the appropriate central geodesic, d — 0. Taylor 
expansion writes 




(19) 
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which we can write as 



w 



^ = ^(;/i^e sin 9 + ^(;/l^e cos 9, 



^(2) 



(3)^ 



(20) 



with 



A* 



e* = (l,0,0) 



(21) 



for j = 2, 3. It will also be useful to define, with j = 2, 3, 



w 



(0)M ^ 

0) dej 



e^ = (l,0,0) '"j' 



(22) 



and 



w 



0) 



e< = (l,0,0)' 



(23) 



so that «;f.) = w'^]^ + wlff. 

Suppose now that the central geodesic of our bundle intersects the worldline of our 
hypothetical source at the point e. The projection of w^^ into the two-dimensional subspace 
orthogonal to both uf^ and to the projection of k^^ into this subspace, i.e. to the direction 
of photon propagation, defines an ellipse. We can determine the characteristics of this 
ellipse explicitly. The relevant projection operator is given by (Kristian and Sachs 1965) 
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HI'. = 5it- 



(24) 



where we denote the ds^ inner product by a dot, e.g. Ue ■ k = u^ga/sk^. We wiU use ± to 
denote the result of acting on a given vector with -ff^i/, e.g. Wj^ = Hf^w'^. 

Extremizing w± ■ w± with respect to 9 we find the major and minor axes of the ellipse 
occur for 



e = -tan"^ 
2 



-2"'^(2) • "'±(3) 



+ 



(25) 



with n = 0,1. In (25) the inner products are evaluated at e. Denoting the argument of the 
inverse tangent by the squared lengths of the two axes are given by, for 9 = (1/2) tan"-^ x 



w± ■ W_i_ = 



2V1 + X' 

- («^±(2) • W_L(2) - W_L(3) • «;_L(3) - 2xW_l(2) • W^±(3)) 



(26) 



and for d = -7r/2 + (1/2) tan-^ x 



W_L-WJ_ = , ^ («^±(2) • «^±(2) + W^±(3) • «^-L(3)) \/l + 

+ («^±(2) • Wj_^2) - Wj_^s) ■ «^±(3) - ^XW±{2) ■ «^±(3)) 



(27) 



which lead to the area of the ellipse, dA, being 
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dA = e^TT 



\ 2 



nl/2 



(28) 



At this point the underlying structure of our calculation is becoming clear. We are 
involved in a realization of standard ideas in the geometry of linear maps of the plane: by 
Jacobi propagation followed by projection we have linearly mapped the two plane spanned 
by w^js^, j = 2,3, to the transverse two-plane at the source. Much of this section is easier 
to understand with this in mind. 

We continue by breaking the projection operator up into zeroth and first order parts 
in the perturbation H^,^ = H^^^^^i, + H^^^^i^. It is then possible to show that to first order 
H^^yw'^ H^aw'^ = H'^^^^^ H^^}w°' . This relationship has a simple geometrical interpreta- 
tion. By a theorem of Sachs (1961), the four- velocity of the emitter has no effect on the 
size and shape of our ellipse so we take the emitter to have v], = 0. In this case H and 
iy(0) project onto two planes inclined with respect to each other by a small angle, the lens 
angle at the emitter. But the projected areas in two such planes differ only by a factor of 
cos a, a the angle of inclination. Since a is first order, the areas are equal to first order. 
Therefore in the expression for the area, (28) above, we are able to replace all i/-projected 
quantities by ff^^^^-projected quantities. 

Next we use the explicit form of the projector to find H^^^^^i^w^ = (0, 0, , w^). Sub- 
stituting back into (28) then yields 



dA = e^TT-l (1 - 20e) Det 

2 



e^7r^(l-2(/)e)DetM^,- 



where the magnification matrix, M^j is given by 
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(29) 



(i)i 



(30) 



with i,j e {2, 3} and 1^ the 2x2 identity matrix. In (29) we have written for r (Ae) and 
7e for 7 (Ae). M*j is a function of the affine parameter and in both the above equations 
is evaluated at A = Ae, corresponding to the point e afong the central perturbed geodesic. 
We may recognize the first of the equalities in (29) as the transformation of area law 
for infinitesimal linear mappings of the plane once we recognize the factors in front of 
the determinant (aside from the factor of tt) as corresponding to the metric factors in the 
induced area form for our image two-plane. Equivalently, this is the induced area two-form 
on the yz— plane acting on H^^^W(^2) and H^^^W(^s^ (again up to the factor of tt). 

The above expression (30) for the magnification matrix is one of the key results of 
this section. We recall from (23) that for i, j = 2, 3, is precisely the variation of the 
transverse separation with the transverse direction cosine. The factor of in (30) turns 
the transverse separation into the transverse angular defiection. Thus, we are begining 
to see the usual structure of the magnification matrix emerge. We notice, however, from 
(14), that while the time variation of the potential does not contribute to the transverse 
deflection of a single ray it will contribute to the gradient in (23) (this will become clear 
below), so that non-static potentials contribute to the magnification matrix. We note 
that (30) may be used for vector and tensor perturbations in addition to the scalar ones 
considered here, after a change in the solution for the separation, hence ^^(j)*, which arises 
from a change in the specific form of the force vector, f^. 

We can use (29) to find d^, the solid angle in the rest frame of the observer defined 
by our bundle. To this end, let Ae = Aq + AA, AA small. We will use (29) to express the 
projected area of our bundle at Ae to order (AA) . Noting that ^sin^^ (Ao — A) = r, we 

have 
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dA (Ao + AA) = Tre^a^ (AA)^ (1 - 2(j)o) 



(31) 



The proper spatial distance in the rest frame of it^ in the metric ds^ corresponding to an 
affine distance Ao; along AL, is given by (Ellis, 1971) AL = |Aa;| (u-k) . bmce 

Ao; = a^AA we can write the projected area of the bundle a unit spatial distance away 
from the observer as 



dQ = dA (AL = 1) 



7re 



[uo ■ ko ) 
= ne' (l + ) (A,)) 



(32) 



where in the last equality we have used u'^ = (l/tto) (l — (po.vl) and = a~^/c^°)^ + 
a~'^k^^^^ = (1 — 20O, —1, 0, 0)). In fact (32) is simply the Lorentz transformation law 
for solid angle to linear order in v^. Were the observer comoving, the solid angle would 
be given by ire^. The expression (32) comes from transforming to a frame in which the 
observer moves with velocity . 

We can use this expression to replace e^7r in (29), yielding 



dA = dQ(l- 2vikf^ (Ao)) ^ (1 - 20e) DetM^j- (Ae) (33) 

This is the desired expression relating the proper area of our emitter and the solid angle 
that it is observed to subtend given its affine distance, Ae (which we will soon trade in for 
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its redshift). The equivalent expression for a congruence of the background spacetime can 
be gained from (33) by taking the perturbed quantities to vanish. This gives 



d^(°) = dl^(°)^ (34) 

where q is some point along the central background geodesic, x^^^^ = (A, r, 0, 0) (see Fig. 
2). 

To proceed we must ensure that the points e and q correspond to the same numerical 
source redshifts in their respective spacetimes. Let the afhne parameter corresponding 
to the point q be given by Xq. We need to impose 1 + z^^^ (Xq) = 1 + z (Ae), where 
2;(°) refers to the redshift to the source in the background along the central background 
geodesic and z refers to the redshift of the source in the perturbed spacetime along the 
central perturbed geodesic. The redshift in the perturbed spacetime is given by 1 + 2; = 
(ue ■ kj (Ae) / (uo ■ (Ao)- 1 + z^^'' niay be found from this expression simply by forcing 
the perturbed quantities to vanish. If we write Xe = Xq + SXq with SXq a first order function 
of Xq, the equal redshift constraint is solved by 



SXq = ^ fefcf ) (A,) - vikl'^ (Xq) +U- aa-'x^'^' + k^'^^] ') . (35) 

Cbq o 

Here a subscript q denotes evaluation at the point q, an overdot denotes an exact conformal 
time derivative along the comoving worldlines, [/]^ = fq ~ fo-, and the separation and its 
affine derivative are solved for along /c^^^)^, which intersects both o and q by construction. 
We have also noted that, to first order, (f)f, = (f)q. We have not used the analogous formula 
for the emitter's peculiar velocity simply to emphasize that the emitter can be moving, 
in principle, under many different types of non-gravitational forces and so the notion of a 
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smoothly varying peculiar velocity field may be inappropriate. 

Recalling the definition of the magnification, we now set dA^^^ in (34) equal to dA 
in (33) and expand Og, Ve, and 7e about their values at q. We also use the equivalence 
of M^j (Ae) and M*j (Ag) to first order. The end result is the following formula for the 
magnification of a source observed at redshift z^^^ (Ag); 



M = 



(DetM^,.) (A,) 



l-2(j)o + 2vikf^ (A,) - 2/c(i)o (A,) 



+ 2 cot«; (Ao - \q) SXq + K sin« (Ao - Xq) e 



(36) 



We emphasize that the source redshift in the actual perturbed spacetime is given 
numerically by z^^^ {^q)- We have chosen to express Ae in terms of Xq rather than the 
other way round because this choice makes it simple to take the physical redshift as the 
independent variable in the magnification formula, (36), rather than Xq. 

An explicit formula for the magnification matrix may be found directly from (30), 
(23) and (14). We need to take a partial derivative of (14) with respect to e*, i = 2, 3, 
while keeping A fixed. The only subtlety arises from the implicit dependence of and its 
spacetime derivatives on e* which arises because of the need to evaluate these quantities 
on a particular unperturbed geodesic. This means 



^ We have, however, considered and Vq to contribute, in a numerical sense, only at 
first order. For instance, the source and emitter four-velocities are properly normalized 
only if Vo(e) are considered to be first order quantities. It is not difficult to relax this 
assumption. 
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d(j) 



= r 



d(f) 

dcj) 



(37) 



In principle 7, and hence gj^J, have a similar implicit dependence on e*, but in fact this 
dependence is vanishing. This is easy to understand once it is realized that varying e* for 
z = 2, 3 with fixed A amounts to variation tangent to the two-sphere at fixed co-ordinate 
radius r(A). Thus r, and so 7 and gjj^J, are stationary. Keeping the above in mind, we find 
an explicit formula for the magnification matrix without restriction to either static or thin 
perturbations (i.e. lenses). With i,j = 2,3, 



M'j = S} + 



sin^ {Xo- Xq 
2 

sin^ {Xo - Xq 



dX sin« (A - Xq)j{X) 



(38) 



X 



5]|^(A)-r(A)7-^(A)^W^^(A)0,,,(A) 



While (38) is computationally quite useful, its relationship to the usual form for the 
magnification matrix is easier to see after rewriting. Introducing the angular diameter 
distance of ds^^\ D (Ai, A2) = a (A2) sin^ (Ai — A2), (see Appendix A), and using the two- 
dimensional projected angle of section 1 above. 



a' = -2^^, 

7 
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(39) 



we can write 



1 .^„.„,,9«- 



y </AO(A,A,)^(A). 



-D(Ao, Ag) ho 

Finally, we note that to zeroth order djde^ = d/dd'^ with ^* the vectorial angle of (1). This 
allows us to make the replacement 



in (40). The result is the most elegant form of our equation. It is also the easiest form in 
which to recover the thin lens limit. For this, suppose the potential is that appropriate to a 
static, geometrically thin lens. Then the first term in (40) vanishes and we can approximate 
the angular diameter factor as constant over the region for which the potential is important. 
Finally, we replace Xq with Ag since, to first order, they are equivalent in M^j. The result 
is 



M-^5;.-£i^A/'^rfAa^ (42) 



We have already seen that the integral of a* over the background path produces the Einstein 
deflection angle. As a result we conclude that in this limit our equation has reproduced 
the usual magnification matrix defined by the ^-gradient of equation (1). 
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5. The Thin Lens in Einstein-de Sitter Space 

As an explicit illustration of the ideas above we will solve for the magnification ma- 
trix appropriate to a point mass perturbation of Einstein-de Sitter spacetime. We will 
find that the usual expression of the magnification matrix is correct provided that the 
impact parameter is much the smallest lengthscale in the problem (save, of course, for 
the Schwarzschild radius of the point mass). Certainly the calculation below is the hard 
way to produce this result. Nevertheless, for the first time the standard result is obtained 
along with correction terms arising from the time variation of the perturbation and the 
difference between the actual path and its piecewise geodesic approximation. 

Our starting point is the potential approximation to general relativity (Martinez- 
Gonzalez, Sanz, and Silk 1990), which writes the perturbation to the fiat FRW metric 
in (5) appropriate to a comoving point particle of mass m located at {xi,yi,zi) on the 
spatial hypersurfaces of constant conformal time as 



It is worth noting that this perturbation is not static. Its time dependence, however, is 
simply that of the background spacetime, that is, it is set by the cosmology. We certainly do 
not expect time variation of this magnitude to affect photons streaming past the lens. We 
will see this prejudice borne out in the calculations below. Nevertheless, this emphasizes 
that more complicated details will emerge in the rigorous picture of lensing than in the 
usual models. 

We will suppose that, other than the point mass, the observer, and the lensed source, 
the spacetime is filled with dust so that we can take a = aorf /r]1 with the subscript o 
denoting evaluation at the observer and Oo constant with the dimension of length (McVittie 
1964). It is useful to keep in mind that with our conventions the only dimensionful numbers 




111 



(43) 



a 
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in this problem are ao and m and both have dimensions of length. 

To avoid unnecessary complications we will take both the source and the emitter to 
be comoving, as would be true, for instance, if both were far from the point mass. We 
suppose the observer to be at the spatial origin of co-ordinates with the source image 
toward the positive x direction. The light from the source does not travel exactly down 
the a;-axis because it is bent by the action of the lens. We could calculate the bending 
by constructing the path using (13) and (14) above. Instead we will use (38) to gain the 
magnification matrix directly. The background path appropriate to the situation is given 
by a;(°)'*(A) = (A, Aq — A, 0,0) with Aq the affine parameter at the observer. We denote 
by Xi that value of the afl&ne parameter at which x^^^f^ intersects the plane x = xi, i.e. 
xi = Xo — Xi- For the usual lensing scenarios, this point of intersection is very nearly the 
point of closest approach of the photon to the lens. 

Given our geometry, the physical linear size of the impact parameter, b, is given by 



It will be convenient to define an unbarred angular-diameter distance symbol, D12, by 
D12 = D{Xi, A2) /a (A2) = (Ai - A2) so that 




(44) 



so that the angle between the lens and the image of the emitter, d, is given by 




(45) 




(46) 



Doi 
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Also, in terms of the vectorial decomposition of 9 introduced in section 1, we have x] = 
6'^ Doll * = 2,3. The physical situation described by the above geometry is illustrated in 
Fig. 3. 

We start the calculation by noting that for an Einstein-de Sitter background we may 
combine the second and third terms in (38) to give 



M%- = 5\ + 



3 "J ■ n 

L>oq -.-o ^^^^ 



The details of our model lens then transform this into 



+ [^2,1 - /o,3 + (Ao + Aq) /i,3 - AoAg/2,3] 5] ^^g^ 

+ [/0,5 - (Ao + Xq) Il,5 + AoAg/2,5] xixj 



aoD, 



oq 



where we have put 



dX- 



n/2 



(49) 



As all of these integrals may be explicitly performed, a general (but complicated and unin- 
structive) expression for the magnification matrix may be written down. For illustration, 
we will consider the situation that the impact parameter is much smaller than all the other 
physical sizes in the problem, so that the integrals simplify to 
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where 



— 

A,5 = 
-^0,5 = 
-^2,3 = 
h,Z = 
Io,3 = 



4 ^ 



2 ^ 

Xie^l, ^xf 



2n2 



Dlq Dol 1 , 

— H — H A 

Xi Xq X, Xn X, 



(50) 



A = 



sinh 



Xi {Xi-xy 

. X9Doi , 



(51) 



With somewhat more algebra it can be shown that /2,i and all the terms proportional 
to A contribute negligibly to the magnification matrix under the assumption of small 
impact parameter. The remaining terms combine to produce 



aoXfDoiDog9^ ' aoXfDoiDoqe^ 



(52) 



which, because q and e are equivalent to the needed order, is also expressible as 
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DolDoe 



( 




) 



(53) 



This agrees exactly with the usual magnification matrix for the point mass lens (Schneider 
et al. 1993, Chapter 2). 

6. Summary 

We have presented formulae (13), (14) for the null geodesies intersecting an observer's 
worldline in an important class of perturbed spacetimes, FRW backgrounds with scalar 
perturbations, in the longitudinal gauge. We have used these equations to obtain a general 
formula (36) for the magnification of ray bundles in these spacetimes. With this, we 
can show for the first time how the usual lens equation (1) and magnification matrix are 
recovered in the curved FRW spacetimes without dividing light paths into near and far lens 
regions. To illustrate our formulae we have calculated the magnification matrix appropriate 
to a point deflector in an Einstein-de Sitter spacetime. We are able to show how the usual 
formula emerges along with (in this case negligible) correction terms. The techniques used 
in this paper are easily applicable to FRW spacetimes with vector or tensor perturbations. 
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Appendix A: The Jacobi Equation of FRW Spacetimes 

We consider the equation of geodesic deviation in the spacetimes ds^^^"^, related by 
conformal transformation to the background FRW spacetimes. The equation of geodesic 
deviation along is 



Here D/dX is the covariant derivative along k^^^i^. We will take to be the wavevector 
of a radial null geodesic parameterized as in section 3 above. Changing variables via 
6x^{X) = P{X,a)^aV°'{X), a an arbitrary affine parameter value along k^^^^^ allows us to 
write (Al) as 



v^'iX) = -KJ^'av''iX). (A2) 



dA2 



We want to solve this equation subject to the initial conditions 5x'^{u) = and dSx^^ / dX{u) = 
5ki^{u). These initial conditions describe some impulse which jolts the geodesic at some 
affine parameter distance u. 

(A2) conveniently decomposes into three separate equations, one for each of the timelike 
component, and transverse and longitudinal projections of the spatial components. To see 
this, multiply (A2) by unity in the form of {S^ — J"pj + J"/?. The result is 
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vi = (A3) 



dA2 



where = (Sp — J*/3^ and Vj_ = J'^ jjV^ . The solutions for the given boundary conditions 
are elementary. Returning to the original variable 5x'^, 



fe°(A) = ^5A:°(«)(A-«) 

7(m) 

5x\{X) = sin«(A - u) 



We can use the solution above to determine the angular diameter distance in the actual 
background spacetime. To do this, let k^^^'^{\) = (1, —7, 0, 0) and suppose an impulse at u 
given by 5k'^{u) = 5k'^{u) = (0, 0, —76, 0), with e infinitesimal. Let 5x'^{X) be the solution 
to the Jacobi equation for this impulse, 5x^{X) = 7(A) sinfe(A — u){0, 0, e, 0) = (0, 0, re, 0). 
Since 5x^ is a Jacobi vector, x(°)'^(A) = (A,r,0,0) and x(o)^(A) + 6x^{X) are neighboring 
null geodesies to first order in e. The angle that their wavevectors make at u is e. The 

1 /2 

proper linear distance (in ds^^^"^) that they span at affine parameter A is (Sx'^gj^JSx'^^ 
The angular-diameter distance of our FRW background is defined as the ratio of this proper 
linear distance to the subtended angle. 



£)(«,A) = ^(fe^^i°W) 

= a(A) sinfj('u — A) 
31 



1/2 



(A5) 



where the sense of the affine parameter in (A5) is that X < u. 
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Appendix B: Geodesic Deviation in Perturbed Spacetimes 

Consider an arbitrary metric perturbed spacetime with metric 



9i^u = gl^J + V (Bi) 



with h^i, smaU. Let x^^^^{X) be an affinely parametrized geodesic of gj^J. Solutions, 
C'^°^^(A), to the Jacobi equation of the background along x^'*'^(A), 



+ 2T%^k^^)'^^^ + r(3;;/c(0)«/c(«)^C^'^'^ = (B2) 

generate nearby geodesies, x^''^(A), via x^^^^{\) = x^^^^{X) + C'-°-*'^(A). Let x'^{X) and 
x^(A) be geodesies of g^^, generated using the perturbative geodesic expansion from x^2^'^{X) 
and x^^^^{X) respectively. We claim that, to first order, a;^(A) — a;^(A) solves the geodesic 
deviation equation of g^i, along Xj^{X). 

There is certainly nothing surprising in this claim. In fact, the assertion must be true 
if the PGE correctly generates nearby geodesies of gf^u as we claim it does. The explicit 
proof offered here can, thus, be thought of as another check on the method itself. Since 
the proof is nothing more than a tedious application of the usual perturbative techniques 
we present it only in schematic form. 

To see that our assertion is true, we start with the Jacobi equation of g^,^ along a;^(A) 



^ + + ^a,,.klk^AW^ = (B3) 
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By a Taylor expansion, the decomposition of the Christoffel terms into zeroth and first 
order expressions in the perturbation, the ansatz = w^^'^^+w^^^^, and the PGE identity 
— k'^A^ + /c*^^-'^, we can express equation (B3) as two equations holding along x^^'^(A), 
one containing only zeroth order terms and the other containing only first order terms. 
The zeroth order equation is seen to be equivalent to (B2), so that we already know that 
t(;(°)^(A) = x^^^^{\) — x^^'^{\) is a solution. The first order equation may be written 
T {w^^'^^ = for some operator T. The specific form of T will not be necessary for our 
current purposes but it is not hard to obtain. 

Consider now the equation obeyed by x^^^'^{X) = x'^{X) — x^^^'^{X), 



d^'xw . ^roiu. (o^cydx^ 



B , 9p(0)Mi.(0)aO^B_ , -pW/x , (0)a, (0)/3 (1)<7 _ ■pCl)^ , (0)a , (0)/3 (^.. 
To. r ad '^B TT r ^ a.B^a'^B i^B ^B — ^ aB '^B i^B \^^) 



^^2 ~'~ Q-P '^B dX al3,a'^B i^B -^B — ^ a/3 

which holds along x2^^(A) (PB). Using Taylor and k^B^^W = k^A^{>)+w^^^^{>)^ where 
• = d/dX^ we can write (B4) as an equation along x^^^{X). If we subtract from this 
equation the equation for a:^''^(A), which also holds along x^"*^ (and is identical to (B4) 
after the subscript -B's are replaced by A^s) we obtain the result T [x'^q — x'^^^ — and 
our assertion is proved. 
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Figure 1. Gravitational lensing by a single symmetric lens. (3 is the angle at the observer, 
o, between the image of the lens, Z, and the unlensed image of the emitter, e. 6 is the 
angle between the image of the lens and the light ray epo, along which the emitter is 
observed. The point of deflection, p, lies in the lens plane. The intersection of the line 
joining the observer to the point of deflection with the source plane is labeled p' . a, the 
deflection angle, is the angle at p between the image of the emitter, e, and p' . D{o,l) 
is the distance between the observer and the point of deflection, D{o,e) is the distance 
between the observer and the source plane, D{l,e) is the distance between the lens and 
source planes. 

Figure 2. The lensing of a congruence of null geodesies. The actual congruence, with 
central ray x^, joins the observer at o with the emitter at e and has area dA at the emitter. 
It is constructed from a congruence of the background, with central geodesic x^^^^^, reaching 
between the observer and a point q. The redshift in the background between the observer 
and q is equal to the redshift in the perturbed spacetime between the observer and the 
emitter. The background congruence has area dA^^^ at q. 

Figure 3. The geometry of Section 4. 9 is the angle at the observer between the lens and 
the emitter. A; is the affine parameter value at which the background geodesic intersects 
the lens plane. Ae is the affine parameter value at which the background geodesic intersects 
the source plane. The observer is located at affine parameter value Aq. Dqi is the angular 
diameter distance in the background between the observer and the lens plane. 
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